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" . Abstract 

fSJ ■ Discussions about Bose-Einstein correlations between decay products of coproduced 

W-bosons again raise the question about the behaviour of correlations if several 
strings are produced. This is studied by the multiplicity dependence of correlation 
functions of particle pairs with like-sign and opposite-sign charge in pp reactions at 

O ! ^ = 630 GeV. 

X 

a: 



1 Introduction 



Recently, there has been much discussion regarding the possibility that Bose- 
Einstein correlations and other interconnection effects between decay products 
of different strings could affect the measurement of the W mass [1]. Since mea- 
surements are hampered by low statistics and experimental difficulties [2-4], 
the question arises whether and where effects of the superposition of sev- 
eral strings can be tested independently. Within the Dual Parton Model [5] 
for hadron-hadron reactions, the charged-particle multiplicity N is expected 
to rise with the number of strings. The decrease in the observed correlation 
strength A of the Bose-Einstein effect as function of multiplicity [6,7] may be 
explained in terms of products of different strings, where each string sym- 
metrizes separately [8,9]. To test this idea further, it is desirable to investigate 
quantitatively the multiplicity dependence of like-sign particle correlations. 

Improvements in experimental analysis techniques [10] and larger data sam- 
ples make it possible to repeat and extend Bose-Einstein analyses with an 
advanced strategy. [11] In this Letter, we investigate correlation functions of 
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particle pairs with like-sign (is) and opposite-sign (os) charges at different to- 
tal charged-particle multiplicities with the same model-independent strategy 
and good statistics. The bias introduced by selecting events of a given overall 
multiplicity is eliminated with the use of "internal cumulants" [10]. 



2 Data sample and normalized density correlation functions 

The data sample consists of 1,200,000 non-single-diffractive pp reactions at 
s/s = 630 GeV measured by the UA1 central detector [12]. Only vertex- 
associated charged tracks with transverse momentum pt > 0.15 GeV/c, \rj\ < 
3, good measurement quality and fitted length > 30cm have been used. To 
avoid acceptance problems, we restrict the azimuthal angle to 45° < \<f>\ < 135° 
("good azimuth"). Since, however, the multiplicity for the entire azimuthal 
range is the physically relevant quantity, we select events according to their 
uncorrected all-azimuth charged-particle multiplicity N. The corrected multi- 
plicity density is then estimated as twice n, the charged-particle multiplicity 
in good azimuth: (dN c /drj) ~ 2(dn/drf). 

All quantities measured are defined in the notation of correlation integrals [13], 



with pj the three-momenta, Pi the corresponding four-momenta, and q = 



\/—(pi — Pt) 2 - As usual, pi denotes the joint particle density of order i which 
in integrated form is the appropriate factorial moment. The integration region 
Q is identical with our experimental cuts as specified above and specifically 
refers to the good-azimuth region. All particles have been assumed to be pions. 

In pp reactions and in full phase space, the number of positive and negative 
particles are equal, as are the corresponding one-and two-particle densities pi, 
Pi- Given the charged particle density pi(p) = Pi~(p)+Pi~(p)i we hence assume 
that pt(p) = Pi~(p) — |pi(p) an d therefore also pt®Pi{Q) — \pi®Pi{Q) 
etc. The like-sign and opposite-sign normalised two-particle densities become, 
respectively, 
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Fig. 1 shows the normalized density correlation functions r 2 for pairs of like- 
sign (£s) charge and for opposite-sign (os) charge separately. Restricting the 
total uncorrected charged-particle multiplicity N in < 3 to the windows 
5 < iV < 9 (Fig. la) and 28 < iV < 35 (Fig. lb), one obtains for the 
corrected particle density in the central rapidity region dN c /dr] = 1.22 ± 0.09 
and dN c /dri = 5.32 ± 0.17 respectively. Both the like-sign and opposite-sign 
correlation densities show a strong dependence on multiplicity.^ 



3 Cumulants for fixed multiplicity and multiplicity ranges 

While the multiplicity dependence of r 2 is of much interest, a quantitative 
analysis must both remove combinatorial background by calculating the cumu- 
lants and also correct for the bias introduced by working at fixed multiplicity. 
This bias arises because at fixed total multiplicity N, standard second-order 
factorial cumulants, defined by 

«2(pi, P2 \N) = p 2 ( Pl , p 2 \N) - Pl ( Pl \N) Pl ( Pl \N) , (2) 

are nonzero even when particles are completely uncorrelated; for example, 
for purely uncorrelated multinomially-distributed events, K™ ult (pi, p 2 \N) = 
— (l/iV)pi(pi \N) pi(p 2 \N) 7^ . Because these correlations result solely from 
the restriction of the sample to a fixed multiplicity, they are termed "external" 
and should be removed. The "internal cumulants" of Ref. [10] 

4( Pl , P2 \N) EE K 2 ( Pl , P2 \N) - KT h (Pl, P2 \N) 

N(N-l) 

= p 2 ( Pl , P2 \N) - { — ' Pl ( Pl \N) Pl ( P2 \N) (3) 

correct this bias exactly: they are zero whenever the iV particles behave multi- 
nomially. 

The following arguments lead to a unique choice of normalization p 2 orm . The 
correctly normalized internal cumulants i^ 2 f ( Pl ,p 2 \N) = K 2 ( P i,p 2 |iV)/p2° rm 
should be independent of multiplicity iV whenever pi(p \N) and p 2 ( P i, p 2 \N) 
depend on N only globally i.e. have the same shapes (in terms of the momenta) 
for different N. Such "shape constancy", 



The usual Bose-Einstein analysis assumes that r 2 tends to a constant for large 
Q, ie. ^ 2 s be(Q > 1) = constant. It should be clear from Figure 1 that no such 
constancy exists for limited-multiplicity windows. 
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N' 

Pi(p\N') = — Pl (p\N), (4) 
p 2 ( Pl , p 2 \N') = N ^ N _^ p 2 (pi, P2 \N) , (5) 

coupled to the requirement that 

K 2 I (p 1 ,p 2 \N') = K 2 I (p 1 ,p 2 \N), (6) 

fixes the appropriate normalisation to be 

P? rm (N) = Pi(Px\N) Pl (p,\N) , (7) 

the quantity to which p 2 defaults when p! and p 2 become statistically in- 
dependent. The correctly normalized internal cumulant at fixed multiplicity 
consequently reads 



JV(JV-l) p,(pi| J V)pi(p 2 |A r ) 

Extending these arguments to multiplicity ranges N e [A, B], specifying them 
for different charge combinations and adopting the correlation integral Eq. 
(1), we arrive at our measurement prescription for the normalized internal 
cumulants for is and os pairs, 



K 2 H \Q\AB) = n \ r e 2 s (Q\AB) - 1 , (9) 

k; Io \Q\AB) = r° s (Q\AB) - 1 , (10) 

n + n_ 



where fi£ S (= n + = n_) , n + n_ and n(n— l) ls are the mean numbers of positive 

or negative particles, os pairs and ++ (or ) pairs respectively in the whole 

interval Q and in the multiplicity range [A, E\. Eqs. (9) and (10) are obtained 
[14] by assuming shape constancy for all averaging procedures^ 

Since all quantities shown here and below are to be understood as mean values 
in multiplicity ranges, we henceforth (and in Fig. 1) omit the bar on the 
symbols. Fig. 2 shows both like- and opposite-sign internal cumulants for two 
selections of dN c /drj. Three features are apparent: 



2 The renormalization factors in front of r 2 were used in Ref [15] and elsewhere. 
Here, they are derived from Eqs. (4), (5) and the requirement (6). 
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I. The importance of changing from r 2 to K\ lies in the fact that the lat- 
ter demarcate clearly the "line of no correlation" which, due to the fixed- 
multiplicity conditioning, is not equal to unity for r 2 in Fig. la. 
II. Internal cumulants integrate to zero over the entire phase space; hence the 
positive part of K 2 at small Q is compensated by a negative part at larger 
Q. Physically, this means that particles like to cluster, so that there is a 
surfeit of pairs at small Q and a dearth of pairs at large Q compared to the 
uncorrelated case. 

III. The dependencies of the like- and opposite-sign cumulants on multiplicity 
are rather similar in that both decrease markedly with dN c /drj. This is 
discussed more fully below. 



4 Multiplicity dependence of normalized cumulants 



The behaviour of K 2 and K 2 os as a function of dN c /drj suggests that both 
could have approximately the same functional dependence C(dN c / drj) on mul- 
tiplicity density. Under this hypothesis, where both depend on dN c /drj with 
the same functional formQ 



K I 2 £s (Q\N c )=Y £s (Q)C(N c ,Q), (11) 
K 2 Ios (Q\N c )=Y os (Q) C(N C ,Q), 

the quotient of the cumulants should be independent of multiplicity, 

Ki ls {Q\N c ) Y es (Q) 



K} OB (Q\N c ) Y°*{Q) 



(constant in iV c ) . (12) 



Figure 3a shows that (12) holds approximately. (In the region where both 
cumulants are near zero, no meaningful quotients can be formed.) 

Having shown that like- and unlike-sign internal cumulants behave approxi- 
mately in the same way as functions of dN c / drj, we now look for an appropriate 
functional form for this dependence. A first hypothesis is that K 2 depends in- 
versely on N c , 

K 2 Ia (Q\N c ) = Y a (Q)C(N Cl Q)=Y a (Q)N; 1 a = £s,os. (13) 
This can be motivated theoretically by 



For simplicity we write N c instead of dN c /drj here and below. 
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a) Resonances: If the unnormalized cumulants k\ os and n 2 is were wholly the 
result of resonance decays and if the number of resonances were proportional 
to the multiplicity N c , then k 2 oc N c . Assuming shape constancy pi(p \N C ) oc 
N c pi(p) gives p\®p\ oc N%, and hence after normalization, the resonance- 
inspired guess is 

/res i 

Pi® Pi N c 

b) Independent superposition in momentum space of v equal strings would 
also lead to K\ = vK 2 string / (pi®p\) oc (1/V) oc (1/N C ). Deviations can occur 
if the strings are unequal or if there is no strong proportionality between v 
and N c . 

Eq. (13) implies that K\ a (Q \N x )jK{ a {Q \N 2 ) = (N 2 /N 1 ) = (constant in Q) 
for two multiplicities N\ and N 2 for the same a (= is or os). In Fig. 3b, we show 
the quotient of cumulants for two multiplicities. Surprisingly, we find not one 
but two regions of approximate constancy in Q, one for small Q < 0.4 GeV, 
one for large Q > 2 GeV, where only the latter corresponds to the value 
(N 2 /Ni) expected from (13), shown as the dotted line. At small Q, one must 
clearly look for other functional forms for C(N C ,Q). Some phenomenological 
guesses are as follows. 

c) In the Quantum Statistical approach of Bose-Einstein correlations [16], 
no multiplicity dependence is expected with our renormalization of r 2 in 
Eq. (9) . 

d) A mixture of processes a) and c) could result in a dependence 

K^(Q\N c )^a(Q) + b -^-. (14) 

However, no comparable picture is available for the unlike-sign case. 

In order to test the above ideas, we plot the cumulants against (dN c / drf) -1 
as follows. To avoid local statistical fluctuations, the normalized cumulants 
K 2 (Q) and K\ os (Q) are fitted with suitable functions in restricted Q-ranges 
for each dN c /dr] (not shown). The best-fit values at small Q (0.1 GeV/c) and 
at large Q (7 GeV/c) are plotted in Figs. 4a and 4b respectively. The Kl ls {Q) 
have also been fitted to an exponential parametrization for Q < 1 GeV/c , 

Kl ls (Q) =a + \e- RQ , (15) 

since the reported increase [6,7] of the radius R in case of is (Bose-Einstein) 
functions could cause part of the decrease with dN c /drj. The A values ob- 
tained are thus corrected for effects of varying radii, but still indicate a pro- 
nounced multiplicity dependence (crosses in Fig. 4a) similar to the model- 
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independent cumulants (filled circles)^} Fig. 4 shows that, as in Fig. 3b, the 
(l/iV c )-dependence is satisfied only for large but not for small Q. The a + b/N c 
dependence in Fig. 4a (solid line) provides a possible, but hardly unique, ex- 
planation. 

The important region around Q — 1 GeV, where the phase space contributes 
maximally [19] is unfortunately difficult to investigate, because there the K\ os 
are decreasing rapidly with increasing Q while the K^ ts are already small 
(Fig. 2). A 1 /^-dependence (due to resonances such as p°) in this dominant 
region around 1 GeV/c could presumably cause the large-Q region to follow 
suit via missing pairs. 



5 Summary 

The multiplicity dependence of like-sign and opposite-sign two-body correla- 
tion functions have been studied with the same model-independent strategy. 
The bias introduced by selecting events of a given overall multiplicity is elim- 
inated by measuring internal cumulants. We observe that 

- the like-sign and opposite-sign cumulants have very similar multiplicity de- 
pendence, 

- there exist two regions, one at small Q, where the multiplicity dependence 
of both is weaker than 1/N C , and one at large Q ( > 2 GeV/c), where the 
cumulants are negative and follow roughly an 1/N C law, 

- a third region arround Q = 1 GeV/c shows small and rapidly changing 
cumulants. 

The decrease of is functions at small Q with increasing multiplicity favours 
an interpretation in terms of a suppression of Bose-Einstein correlations be- 
tween products of different strings. In the Dual Parton Model approach [5], 
multiparton collisions corresponding to multipomeron exchange are expected 
to contribute to the inelastic cross section. The observed energy dependence of 
multiplicity distributions supports this view [20] . Up to 2-3 pomeron exchanges 
might occur at our highest multiplicities. This could explain quantitatively the 
corresponding suppression of is (Bose-Einstein) functions in Fig. 4a. But this 
interpretation is not unique. Adopting e.g. the "core-halo picture" [21], one 
could explain the observations by assuming that long-lived resonances (such 

4 The values of R from the fit to Eq. (15) increase by about 30% over the range 
of dN c /dr] considered. A stronger dependence of R on multiplicity and strongly 
decreasing A values have been observed when fitting Gaussian functions [6] and 
extending the fit range to 2 GeV/c. This indicates the dependence of the results on 
the choice of fit functions and regions. Here, we are emphasizing the region of small 
Q where Gaussian fits fail completely [17,18]. 
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as the uo and 77) are produced more frequently at larger N c thereby increase 
the relative strength of the halo [22] and hence do not contribute to Bose- 
Einstein correlations. Resonance decays (Regge terms) should arguably also 
contribute to £s-functions [23] . One could therefore try to explain alternatively 
the decrease of is functions with iV c by a mixture of Bose-Einstein correlations 
(assumed to be constant in N c ) with resonance production. 

All this leaves unexplained, however, the similar behaviour of os and is func- 
tions. Theoretical work [5,8,9,24] is challenged by the above experimental re- 
sults. Together with the results of other reactions, they represent a new piece 
of information [25] in the colourful puzzle of multiparticle production. 
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Fig. 1. Normalised moments r| s versus Q for like-sign pairs (filled circles) and 
rf versus Q for opposite-sign pairs (open cirles) at a) dN c /dr] = 1.2 and b) 
dN c /dr] = 5.3. 
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Fig. 2. Internal cumulants Eqs. (9)— (10) for is pairs (filled circles) and os pairs 
(open circles) for multiplicity densities a) dN c /dr] = 1.2 and b) dN c /dr} = 5.3. 
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Fig. 3. a) The ratio Eq. (12), shown for two selections of dN c jdr\. b) Ratio of 
two like-sign and two opposite-sign cumulants at different multiplicities dN c /dr). 
The region around Q = 1 GeV has been omitted because both numerators and 
denominators are near zero. 




Fig. 4. a) Multiplicity dependence of K^ s (filled circles) and K^ 08 (open circles), 
both at Q = 0.1 GeV, plus A values (crosses) of exponential fit (15). b) as in a) but 
for the large Q = 7 GeV-region. Note from Fig. 2 that the cumulants are negative 
at large Q. For better comparison of the respective dependencies on dN c /drj, the 
absolute values have been scaled by constant factors. Solid lines are best fits to the 
A's using Eq. (15). 
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